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Hybrid inflation is a two-field model where inflation ends due to an instability. In the neighborhood 
of the instability point, the potential is very flat and the quantum fluctuations dominate over 
the classical motion of the inflaton and waterfall fields. In this article, we study this regime in 
the framework of stochastic infiation. We numerically solve the two coupled Langevin equations 
controlling the evolution of the fields and compute the probability distributions of the total number 
of e-folds and of the infiation exit point. Then, we discuss the physical consequences of our results, 
in particular the question of how the quantum diffusion can affect the observable predictions of 
hybrid infiation. 
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I. INTRODUCTION 

Inflation is the leading scenario among the models at- 
tempting to describe the physical conditions that pre- 
vailed in the very early Universe. It consists in a phase 
of accelerated expansion which naturally solves the prob- 
lems of the hot big bang theory [ll-Q (for reviews, see 
Refs. [E-dl)- In addition, it predicts an almost scale in- 
variant power spectrum for the primordial cosmological 
fluctuations, the tiny deviations from scale invariance be- 
ing related to the microphysics of inflation (9l-[l3| . As is 
well known, this prediction turns out to be fully consis- 
tent with different types of astrophysical observationsm 
among which is the measurement of the cosmic mi- 
crowave background radiation (CMBR) anisotropies [isl - 

Inflation is usually driven by one or many scalar fields. 
In the context of general relativity, this represents the 
simplest mechanism to obtain the negative pressure nec- 
essary to produce an accelerated expansion. However, 
the physical nature of those scalar fields is presently 
unknown, and many different inflationary models have 
been suggested. The reason for such a situation origi- 
nates from the fact that inflation is a high energy phe- 
nomenon. Indeed, the energy scale of inflation is some- 
where between the TeV scale and the grand unifled the- 
ory scale [15]. At those scales, particle physics remains 
elusive, and for a consequence, there is presently a large 
variety of different inflationary scenarios. 

However, given the extensions of the standard model of 
particle physics, notably those based on supersymmetry 
and supergravity, it is clear that some models appear to 
be more motivated and more generic than others. In par- 
ticular, this is the case of hybrid inflation ^OjHH, which 
can be realized in various ways in the context of super- 
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symmetry, see for instance the scenarios named i^-term 
inflation and Z?-term inflation (among others) (22l - [25| . 
Hybrid inflation is a two-field model such that inflation 
occurs along a valley in the field space and ends by tachy- 
onic instability along the so-called waterfall field direc- 
tion. Hybrid infiation is known to lead to a blue spec- 
trum for the fiuctuations, a prediction which appears to 
be disfavored by the most recent observations [15] . How- 
ever, it was shown recently p6l - l28l | that, in some regions 
of the parameter space, a significant number of e-folds 
can occur in the waterfall regime. In this case, it was 
also demonstrated that the spectral index becomes red, 
which therefore implies that the model is in fact totally 
compatible with the data [26l, • 

In the context of inflation, another interesting ques- 
tion is the role played by the quantum corrections [29l - 
[39| . Various works have shown that they can have a 
crucial impact on the inflationary dynamics. This is, for 
instance, the case for large field infiation if one starts in- 
flation high enough in the potential. In this case, the 
quantum kicks undergone by the fleld can be so impor- 
tant that the field climbs its potential instead of rolling 
down it as should be the case according to the classical 
equations of motion. In such a situation, it is likely that 
one enters into a regime of eternal inflation [sil, Ho, |4l| . 

Hybrid infiation is also a model where one expects the 
quantum corrections to be very important. It should be 
the case high in the inflationary valley but also around 
the critical point where the tachyonic instability is trig- 
gered d^, [13, ■ The goal of this article is to investigate 
this last question in detail. In particular, we are inter- 
ested in whether the quantum effects can significantly 
modify the classical dynamics and affect the observa- 
tional predictions of the model. 

In order to carry out our s tudy , we use the stochastic 
inflation formalism [sol. l32| - [39l l43l . It consists in mod- 
eling the quantum effects by a stochastic white noise. As 
a consequence, the equation describing the motion of the 
flelds becomes a Langevin equation. As mentioned previ- 
ously, hybrid inflation is a genuine two-fleld model, which 



implies that one has to deal with two coupled Langevin 
equations. Solving this system is a very complicated task, 
even if a perturbative expansion is used, as usually done 
in the context of single field inflation. This is the rea- 
son why, in this article, we use a numerical approach. 
This allows us to compute various interesting quantities 
such as the probability density function for the number 
of e-folds or for the location in field space of the end of 
inflation. 

This article is organized as follows. In the next section. 
Sec. ini we review in some detail the classical behavior of 
the inflaton and waterfall fields. This allows us to clearly 
identify the region in the parameter space where a sig- 
nificant number of e-folds can occur during the waterfall 
regime. This also permits a comparison between the clas- 
sical and stochastic dynamics. In Sec. lIIII we numerically 
solve the two coupled Langevin equations that control 
the behavior of the two fields. We then use this result 
to compute various probability density functions, in par- 
ticular, that of the number of e-folds and of the inflation 
exit point. Finally, in Sec. IIVI we summarize our main 
results and present our conclusions. 



II. CLASSICAL REGIMES 

There exist many ways to realize hybrid inflation. In 
this article, for simplicity, we focus on the first version 
studied in Ref. ^2(J]; see also Ref. [1^. In this case, the 
potential in the field space (0, ip) , where (j) is the infla- 
ton and -ip the waterfall field, is given by the following 
expression: 



1 - 



M2 



02 ^^02^2 
^2 



02 M2 



(1) 



It contains four parameters (of dimension one). A, M, /i, 
and (pc- The scale A is fixed by the Cosmic Background 
Explorer (COBE) normalization (the other parameters 
being fixed) . The true minimums of the potential are lo- 
cated at = and ip = ±M, while the instability point 
is given by = 0c, "0 = 0. Along the inflationary valley, 
tjj = 0, the potential reduces to A'' [l -I- (0//i)2] which 
shows that, in this regime, inflation cannot end by viola- 
tion of the slow-roll conditions. The full hybrid inflation 
potential is shown in Fig. [1] where the inflationary valley 
is clearly visible. 

In this section we study the classical behaviors of the 
inflaton and waterfall fields. The slow-roll equations con- 
trolling the evolution of the fields can be expressed as 
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where H = a/a is the Hubble parameter, a{t) being 
the Friedman-Lemaitre-Robertson- Walker (FLRW) scale 
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FIG. 1: Potential of hybrid inflation in the (4>,'4>) plane. The 
values of the parameters are /i = 3190.4 Mpi, M = 0c = 
0.1503 Mpi, with Mpi being the reduced Planck mass. 



factor and a dot denoting a derivative with respect to 
cosmic time. The quantity N is the number of e-folds, 
N = ln(a/ai), where a{ is the scale factor at the begin- 
ning of inflation. 

In order to study the classical dynamics, it is interest- 
ing to calculate the slow-roll parameters. The hierarchy 
defined from the potential |46l - l49{ is given by the follow- 
ing expressions: 
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On the other hand, the hierarchy defined from the Hubble 
parameter, the so-called Hubble flow parameters [50..,51il. 
can be expressed as 
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where eg = H-J H{N). The above expression implies that 
having inflation is strictly equivalent to ei < 1, where 
€i = —H/H^. Obviously, the two hierarchies are related 
by simple expressions. In particular, the first horizon 
flow parameter is 



ei — e<^ + ev 
where and have been defined before. 
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FIG. 2: Exact numerical solution for the inflaton (left panel, solid blue line) and waterfall (right panel, solid blue line) fields 
in the inflationary valley . The red dotted lines represent the analytical solution and it is obvious that the approximation is 
very good. The damped oscillations of the waterfall field along the inflationary valley are clearly visible. The WKB analytical 
formula given by Eq. (|17|) is a very good approximation to the exact numerical solution. 



Having specified the notation, we now turn to the 
choice of the free parameters controlling the shape of the 
inflationary potential. 



A. Physical Priors 

It is usually assumed that hybrid inflation occurs in the 
vacuum dominated regime, for which cj) ^ fi and ^ <^ M . 
In this paper we also assume that this is the case. In 
any case, otherwise, hybrid inflation in the valley would 
be equivalent to a large field model, which is not the 
regime of interest here. For simplicity, in order to reduce 
the number of free parameters, and as also motivated by 
the supersymmetric version of the model, we take (j)c — 
M. Notice that this assumption does not imply a loss of 
generality, as we could easily relax it without drastically 
modifying the results obtained in this paper. Finally, 
in order for inflation to proceed for small values of the 
fields (compared to the Planck mass), one can consider 
that 4>c, M <^ Afpi, Afpi being the reduced Planck mass. 

In the valley and in the ^ 1 limit, the slow-roll 
parameters ei and 62 read 
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Therefore, for the slow-roll approximation to be satisfied, 
these two parameters have to be much smaller than 1, 
which implies that n 3> Mpi. 

In the next subsection we study the behavior of the two 
fields during the first phase of evolution, namely when the 
infiaton field is moving along the valley. 



B. The Inflationary Valley 

The question of the initial conditions in hybrid infla- 
tion is a very interesting and non-trivial question. It has 
been studied in detail in Refs. [26l. [s^ - ls^ l . Here, we sim- 
ply argue that starting in the valley can be reasonably 
justified even if more complicated regimes can be found; 
see Ref. [HI]. Indeed, if inflation starts beyond the critical 
line 4> = (j)c, the system very quickly reaches the region 
where ip/M ^ 1. In this regime, where the infiaton field 
is driving inflation, the slow-roll equation of motion for 
(j) can be integrated, leading to 
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where (j)-^ denotes the initial value of the infiaton field. 
This relation can be inverted, and one obtains |15l] 
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where Wq denotes the 0-branch of the Lambert function. 
In the <C 1 limit, this formula simply reads 
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This last expression is compared with an exact numerical 
integration of the full equations of motion in Fig. [2] (left 
panel), where it is shown that this is indeed an excellent 
approximation. Moreover, this allows us to calculate the 
number of e-folds "generated" in the valley, which reads 
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Clearly, this number is large because /i ^ Afp 
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FIG. 3; Exact numerical solution for the waterfall field 
(solid blue line) in the inflationary valley after the oscillatory 
regime. The red dotted line represents the analytical solution, 
and it is obvious that the approximation is very good. 



Let us now study the behavior of the waterfall field 
in the vicinity of the valley, when tp/M <^ 1. Since ip 
undergoes damped oscillations in this regime, it is clear 
that the slow-roll approximation cannot be used. On 
the other hand, since "tp oscillates much faster than p 
moves, the Wentzel-Kramers-Brillouin (WKB) approxi- 
mation can be used to describe this regime. The solution 
can be expressed as 



tjj = ipi„ 
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where ipij^ is the initial value of the waterfall field, and 
where / (N) is defined by the following expression (an 
unimportant sign has been ignored): 
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and where Ci and C2 are integration constants. The 
validity of the WKB approximation can be checked 
by estimating the following quantity, (do; / dA'^) /w ~ 
0(l)Mp|//i2 ^ I From the above expression one notices 
that oscillations in the tp direction occur in the regime 



16M2 <^2/^2 



1 



3 M2 l-|-02/;i2 

that is to say, in the region 

<P 3 M2 

— > 1 



> 1 



32 Ml 



(20) 



(21) 



Since M <C Mpi, we see that the field oscillates almost 
all the time before the critical point is met. In fact, it 



turns out that the integral I{N) can be performed. One 
finds 



I{N) ~ -V3. 
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This solution is compared with the exact numerical solu- 
tion in Fig. [5] (right panel). Clearly, the approximation 
is excellent. 

When the condition (|2T|) is not satisfied, I{N) ~ 
±3A^/2 + ■ ■ ■ and the oscillations stop. Since the gra- 
dients become small, this time, one can use the slow-roll 
approximation in order to describe the motion of ip. No- 
tice that, since M is small (in Planck units), the above- 
mentioned regime occurs for a very small range of values 
for (p. However, a large number of e-folds cx fi^ M'^ / M^^ 
can be produced during this phase. The slow-roll equa- 
tion of motion for ip can be straightforwardly integrated 
and gives 
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where p (n) is given by Eq. (IHl) . Since (p > pc the 
valley, ■0 decreases with TV and obviously remains in the 
Ip ^ M region. If one uses the fact that (p ^ then the 
integral in the above formula can be performed exactly. 
Upon using Eq. ([T5|) . one obtains 



tp = ipinexp 
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where, in the last equation, we have used the fact that 
fj, ^ Mp|. In that case, one concludes that is exponen- 
tially damped after the oscillations have stopped and be- 
fore the critical point is reached. To our knowledge, this 
regime was not considered before. The above expression 
is compared with an exact numerical integration of the 
full equations of motion in Fig. [3l As one can notice, 
the agreement between the exact numerical solution and 
the analytical approximated expression is excellent. The 
previous formula also allows us to calculate the classical 
value of ■0 when the system reaches the critical point. It 
is given by 



ipc = ipin exp 
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In practice, this value is always extremely small, thanks 
to the fact that fi Mp^ and, as was noticed in Refs. [13, 
|4^ . the quantum fluctuations of tp can be much larger 



5 



than its classical value. We come back to this point in 
detail in the next section. However, before addressing 
this issue, in the next subsection, we describe the classical 
motion of the two fields during the waterfall stage. 



Then, using the slow-roll trajectory one easily finds that 



In- 



0c 8^2^,2 



1 



f 



(32) 



C. The Waterfall Regime 

The waterfall regime has rece ntly been studied by var- 
ious authors; see e.g. Refs. [27l. [SSt^SSl . Here, we mainly 
follow the terminology used in Ref . [Sg] . We assume that 
slow roll is valid initially, at the critical point. We first 
study the so-called "phase 0" [1^ . It consists in neglect- 
ing the last term in the inflaton slow-roll equation ([2]) 
and the first one in the right-hand side of the waterfall 
equation ([3]) (on the grounds that, initially, = 0c). No- 
tice that, in this case, one could also solve the full infla- 
ton equation, keeping the second time derivative, since 
in this approximation it becomes linear. On the other 
hand, the waterfall equation is nonlinear. In this sense, 
we do not start from a linear situation. It is easy to find 
the (slow-roll) solutions, and they read 
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where denotes the number of e-folds at the critical 
point, i.e. at the onset of the waterfall phase. In field 
space the trajectory reads 
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Of course, instead of expressing </> in terms of ip, one can 
also express in terms of the inflaton field. In this case 
one obtains 
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These expressions are fully consistent with Ref. [58|. 

The next question is when phase stops. By definition, 
upon using Eqs. ^ and ([3]), it occurs when = 0i and 
■0 = "01 such that 
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Indeed, among the two conditions that we have required 
in order to derive the solutions ([25)1 and (|27p . this one is 
the first to be violated since ip{N) decreases during phase 
0. This condition can also be written as 
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If we are in the regime where S^'^tp'^/M^ <^ 1, then one 
has 
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From these expressions one can easily estimate the num- 
ber of e-folds in phase 0. One obtains 
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where Ni denotes the number of e-folds at the end of 
phase 0. We see that the above quantity [as well as the 
parameter used in the expansion that leads to Eqs. p4|) 
and ([55])] is controlled by fi/AIpj, which is large, and by 
-0c/M which is small. Therefore, the smallness of this 
parameter is a priori not obvious. The two situations, 
where it is small or large, have been studied in Ref. [ssj . 
However, in practice, ipc/M is so small that the parame- 
ter mentioned previously is always small. In this case, we 
conclude that phase is unimportant since it lasts a neg- 
ligible number of e-folds. As a consequence, the values of 
and "0 remain almost unchanged during that phase. 

We now proceed with phase 1. By definition, the sec- 
ond term on the right-hand side of the waterfall equa- 
tion ^ can be neglected. This means that this equa- 
tion, as was already the case for the inflaton equation of 
motion (which remains unaffected during phase I), be- 
comes linear. For this reason, sometimes, this phase is 
also called the "linear phase." During this phase, the so- 
lution for the inflaton field is unchanged but, of course, 
one now has to solve the new approximated equation for 
the waterfall field . The solution can be easily calculated 
and reads 
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Then, one can Taylor expand the exponential functions 
since we are in the regime where /^/Mp, ^ 1. This gives 
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FIG. 4: Exact numerical solution for the inflaton (left panel, solid blue line) and waterfall (right panel, solid blue line) fields. 
The red dotted lines represents the slow-roll analytical solution during phase 1. 



This solution together with the solution for the inflaton is 
represented in Fig. |4l This plot confirms that the previ- 
ous approximated solutions match the exact ones with a 
very good precision. Finally, in field space, the trajectory 
now reads 



ip = ipi exp 
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or, equivalently, 
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Phase 1 stops when, on the right-hand side of the slow- 
roll inflaton equation of motion ([2]) for (p, the last term 
becomes 1. This occurs for ijj = ^p2, where 
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the last approximated relation being obtained under the 
assumption that phase can be neglected and, as a con- 
sequence, that 01 ~ 0c and ipi — ipc- It is also important 
to realize that the terms 1 — 02/02 tp'^/M'^ are equal 
at the onset of phase 1 and then both increase. It is 
therefore necessary to check that, at the end of phase 
1, the term 1 — 0^/0c still dominates over ip'^/M^. In 
other words, one has to verify that ■02/^2 has increased 
less rapidly than 1 — 0^/0c- Using the solution for the 
waterfall, one has 
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where N2 denotes the number of e-folds at the end of 
phase 1 or, equivalently, at onset of the phase 2. Upon 
using this formula, this leads to 
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an expression that should be compared with 
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We see that the condition 0|/0c — 1 ^ ij^l/M'^ is a priori 
not obvious. However, in the case under scrutiny in this 
article, one chooses 0c and M to be roughly of the same 
order of magnitude and /i ^ Mp,. As a consequence, 
the condition is satisfied since 02 /0c — 1 scales as the 
inverse of n (neglecting the influence of the logarithm) 
while ?/;|/M2 scales as the inverse of . However, it is 
also clear that one could easily design a situation where 
this is not true. Here, we restrict ourselves to situations 
where this does not happen. 

Finally, let us express the number of e-folds produced 
during phase 1. It is given by 



^2 - iVc 



2^/2M^ 



In 



1/2 f 4>cM 



(47) 



Upon using Eq. ([25|) . one could also replace ipc by its 
expression in the above equation to obtain a formula de- 
pending on the initial conditions. We see that the number 
of e-folds during phase 1 is essentially controlled by the 
ratio /iM/Mpj. This conclusion is in agreement with the 
results of Ref. [11]. As a consequence, for /xM/Mp^ > 1, 
-^2 — ^c can be large. Hence, we conclude that the num- 
ber of e-folds during the waterfall phase can indeed be 
much greater than the 60 required for inflation to be suc- 
cessful as was noticed in Ref. [l^l- The previous con- 
siderations allow us to identify where, in the parameter 
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FIG. 5: Classical background evolution of the inflaton and 
waterfall fields starting from the critical point, <j) = <j)c and 
■4) = fO"^^Mpi. The solid blue curve represents the exact 
(i.e. numerical) trajectory while the dotted red curve is the 
slow-roll approximation during phase f and the dotted green 
curve is the slow-roll approximation during phase 2. The 
inset shows the overall evolution of the two fields with, in 
particular, the oscillations around the true minimum of the 
potential at the end of inflation. 



position where we can turn to the main topic of this ar- 
ticle, namely, studying the role played by the quantum 
effects. This is the goal of the next section. 



III. STOCHASTIC EFFECTS 

In this article, we use the stochastic inflation formalism 
to study the quantum effects. In this formalism, an ef- 
fective Langevin equation can be derived for the "coarse- 
grained" field, i.e. the original field averaged over a phys- 
ical volume the size of which is typically larger than the 
Hubble radius . Applied to the case of hybrid infla- 
tion, one obtains two coupled Langevin equations for the 
inflaton and the waterfall fields, respectively. They read 
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where and arc two uncorrelated white Gaussian 
noises with 0-mean and 1-variance. Notice that the time 
variable used is the number of e-folds. It was argued in 
Refs. [59l - [6l| that this choice is preferred. 



space, this regime occurs. We have studied this classical 
phase of evolution in some detail because this regime is 
of particular interest for the present article. Indeed, in 
the next section, we show that in this case the quantum 
effects play an important role. 

Let us now briefly mention phase 2 (it was studied 
in more detail in Rcf. 58]). This time one needs to keep 
the last term in the inflaton equation of motion ^ . This 
means that the evolution for (p is modified and, as a con- 
sequence, the formula giving 'iJj{N) is no longer valid since 
it made use of the evolution for established before. In 
this regime the Eq. ^ and ^ become fully coupled. 
However, it is still possible to find the slow-roll trajec- 
tory in the field space. One obtains j58i] 



dip 02 _ 02 

which can be easily integrated, and the solution reads 



(48) 



202 In 

02 



(49) 



This expression (green dotted line) is compared to the ex- 
act numerical solution (blue solid line) in Fig. \5\ Clearly, 
the agreement is excellent. During phase 2, inflation 
stops and the system starts oscillating around one of the 
two true minimums of the potential. This is the onset of 
the reheating phase. 

The above considerations complete this section. Hav- 
ing mastered the classical dynamics of the fields in the 
valley and during the waterfall regime, we are now in a 



A. Can the Quantum Effects Be Important? 

Having at our disposal the two Langevin equations pre- 
sented above, the first question is whether the stochas- 
tic noises can really play an important role and, if so, 
where in the field plane. This issue can be addressed 
in the following manner. During a typical time interval 
At = H~^, both stochastic and classical evolutions of the 
fields and tp can be read off directly from Eqs. ([SOt and 
dnH). Roughly speaking, the typical classical change in 
the inflaton value is ~ Mp^{dV/d(j))/V, while the magni- 
tude of the quantum kick is H/{2t:). Therefore, in order 
to assess the relative contribution of the stochastic effects 
over the classical ones, one can study the ratios and 
of these two quantities for each field (in the context 
of a quartic large field model, this is how one can deduce 
that the quantum corrections dominate if the value of the 
field is larger than A~^/^Mpi, where A is the self-coupling 
constant that appears in the potential). This amounts 
to taking A^ ee V^/^ /[2TrV^M^,{dV/d(j})] and a similar 
definition for A.^. In the vacuum dominated regime, the 
two quantities A^ and A.^ read 



1 A2 



A^ = 



47rV^ Mp3 
1 A2MM 



^ 



'M2 



1 



(52) 



(53) 



These quantities are plotted in Fig.lBlin the (0, ip) plane. 
Values such that A > 1 indicate that the quantum effects 
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dominate. Let us now discuss how and A^ behave 
in the field plane. Clearly, A^ is infinite when (jj = 
for any values of ^p. Therefore, the quantum effects are 
dominant along that direction. In the inflationary valley 
'0 = 0, which is perpendicular to the previously men- 
tioned direction, one has 



L valley 



1 

This means that A^ > 1 as long as 

^valley I ^ 



Mp 



< 



Mp 



47rV3 M4 



(54) 



(55) 



For the parameters used in Fig. [6l one has 0™"°y/Afpi ~ 
94, i.e. a value much larger than the upper limit of this 
plot. This means that the quantum effects dominate 
"very high" in the valley and, in particular, around the 
critical point. The previous considerations explain the 
cross-shaped white region centered at the origin observed 
in Fig. [5] In this regime, one expects a faithful descrip- 
tion of the system to be obtained only if the stochastic 
noises for the two coupled fields are taken into account. 
In the following, we study this case, where treating one 
field (for instance, the inflaton) classically and the other 
(the waterfall field) stochastically is a priori not a good 
approximation. 

Of course, these results also depend on the parame- 
ters, in particular, on A. It is interesting to determine 
the value of A = A^ such that (/j^''^^ = 0^- This value 
indicates the limit between the regime where it is manda- 
tory to take into account the noise both in the inflaton 
and waterfall field directions and the regime where the 
waterfall field is still stochastic but where it is sufficient 
to treat the inflaton classically. It is given by 

Al = 4.V3^ 



(56) 



For our fiducial parameters, this leads to A^ ~ 5.7 x 
10~'*Mpi. Obviously, for larger values of /i, A^ is even 
smaller. In Fig. [71 we have represented A^ for the same 
parameters, except that A = 5 x 10~'*Mpi < A^. This 
plot confirms that the region A^ > 1 covers a much 
smaller area which does not encompass the critical point. 
In that case (j> should behave almost classically in the val- 
ley, and we will also investigate this regime in Sec. IIIIDI 
Let us now describe A^. It is infinite for ip = 0, that is 
to say, in the valley. When V' 7^ 0, the quantum effects are 
dominant when (j>^ /(j>l ~ 1 — ijp' /M"^ , i.e. in the direction 
(/) ~ perpendicular to the valley (in the regime where 
<^ 1). This explains the cross-shaped white region, 
this time centered at the critical point; see Fig. [6] This 
time, the previous considerations do not depend on A, 
which means that the noise in the waterfall field should 
always be taken into account in the valley and around the 
critical point. Since this corresponds to a very flat region 
of the potential where most of the e-folds are realized, 
one can already expect the inflationary dynamics to be 
signiflcantly affected by the quantum effects. 



B. Obstacles to a Perturbative Approach 

Having justified that the quantum corrections play a 
crucial role, the next question is how to compute them, 
i.e. how to solve the two Langevin equations. It is clear 
that an exact analytical solution is not available. How- 
ever, as proposed in Refs. [i^li^l, the Langevin equation 
can be solved perturbatively by considering the coarse- 
grained field as a perturbation on top of the classical 
solution. The corresponding formalism in the case of sin- 
gle field inflation was presented in Ref. [i^. However, 
in the present case, we are in a two-field situation, which 
means that both the infiaton and the waterfall fields must 
be expanded according to 



0ci (N) 
V-ci (N) 



S(t>i (N) - 
- (501 (N) 



S(t>2 (N) - 
- 6^J2 (N) 



(57) 
(58) 



where 0ci and V-'ci are the classical values. We see that 
the corrections to the classical solutions are obtained by 
adding successive terms of higher and higher powers in 
the noise. In Ref. [43], general formulas, valid at sec- 
ond order, are provided, leading to a Gaussian proba- 
bility density function for the field. The validity of this 
approach relies on the smallness of the stochastic effects 
compared to the classical ones and, obviously, the expan- 
sion is valid only in a limited regime; see Ref. |44| . Here, 
we have just seen that the quantum effects are dominant 
around the critical point and, therefore, there are already 
reasons to guess that a perturbative approach is not very 
appropriate. 

Moreover, one can see that the perturbative approach 
is technically impossible to carry out in a multiple field 
situation since even the linearized coupled stochastic dif- 
ferential equations cannot be analytically solved. Indeed, 
in the hybrid inflation case, at first order in the noise, 
they can be written as 
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(59) 



(60) 



where is the derivative oi H — W / (Mp,\/3) with 
respect to </> and the other notations used in this equation 
straightforwardly follow. The matrix of this differential 
system does not commute with itself at different times 
N and, as a consequence, one cannot solve the coupled 
perturbative problem in a simple way. 
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FIG. 6: and plotted in the {4>, ip) plane for the parameters n = 3190.4 Mpi, M ^ cj>c ^ 0.1503 A'/pi, A = 0.01418 Mp,. A^ 
and/or A.^ greater than 1 (~ in white on the plot) signal that the quantum effects are dominant. The stochastic effects in the 
^ direction obviously dominate over the classical contributions in the valley and around the critical point, while the stochastic 
effects in the direction dominate in the valley and around the origin. 
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FIG. 7: A^ plotted in the (0, ip) plane for the parameters 
At = 3190.4 A/pi, Af = 0c = 0.1503 Afpi, A = O.OOOSATpi < A^. 
For this value of A, A^ remains small along the valley and 
does not encompass the critical point. It becomes larger than 
1 only in the vicinity of the origin. 



FIG. 8: Numerical predictions for (^/'c)num normalized to 
(V'c)port given by Eq. (|6ip for different values of A normal- 
ized to A^ defined in Eq. (|56p . 



C. Testing the Numerical Approach 

For all these reasons, the only method left seems to be 
a full numerical integration of the stochastic inflationary 
equations. This is the method used in the present ar- 
ticle. Since the differential equations to be solved turn 
out to be stiff most of the time, we use a fourth order 
Rosenbrock method, monitoring a local truncation error 
to adjust step sizes, that we have adapted to take into 
account the presence of an extra random stochastic term. 
When possible, we have also used the Euler-Muruyama 
method in another independent code in order to check 
our numerical results. 



In this section, we describe the tests that we have per- 
formed in order to check that our numerical codes work 
properly. A first verification of the consistency of our nu- 
merical treatment can be obtained in the following man- 
ner. If one considers that the dynamics of (j) remains 
classical in the valley, then following Refs. [Ij, one 
can perturbatively estimate the typical dispersion of the 
waterfall field distribution (i.e. in the ijj direction) at the 
critical point. One obtains 



pert 



327r3/2A/2 



(61) 



Therefore, numerically, in the regime A < A^ (to ensure 
that the inflaton field behaves classically), one should 
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FIG. 9: Stochastic trajectories in field space for different values of A: A = A0/2O (top left panel), A — A^/10 (top right panel), 
A — A^ (bottom left panel), and A = lOA^ (bottom right panel). The solid black line represents the classical trajectory starting 
from the point (j>in ~ 1.0001(^c and tpin = 6.6 x 10~^A/. The blue, cyan, pink, and purple lines represent four different stochastic 
trajectories. The parameters chosen are ji/Mpi = 3190.4 and M/Mn = (pc/Mm ~ 0.1503. The contours A0 = 1 (solid red line) 
and = 1 (solid green line) are also represented. 



recover the same result. As a consequence, it is inter- 
esting to plot the quantity ((V^c)num/(V'c)port)^^^, where 
('0c)num is the dispersion in the waterfall direction (at 
(j) = 4>c) obtained numerically. This quantity as a func- 
tion of A/A^ is represented in Fig. [51 As it is clear from 
this plot, when A < A^, the ratio (('0c)num/(V'c)port)^/^ 
is precisely 1, thus showing that our code correctly repro- 
duces the known analytical result. We also see that when 
A > A^, the perturbative regime breaks down. From 
what we have just discussed, the interpretation of this 
result is clearly that the stochastic effects in the (f) direc- 
tion play a role and kick the system below the critical 
point more rapidly. As a consequence, the distribution 
in ip has much less e-folds to broaden than classically 
predicted, and hence (V'?)num < (■'/'?)port- This behavior 
is similar to what has been found in Ref. [6^ . where it 
has been shown that, in case of a multiple field inflation- 
ary dynamics with one flat direction and several nonflat 
directions, the fluctuations of the nonflat directions can 
be sufficient to block the growth of the root-mean-square 
amplitude along the flat direction. The fact that (V'c)num 



deviates from (V'c)port precisely at A = A^ is another in- 
dication of the consistency of our numerical results. 

Another type of consistency check can also be per- 
formed by investigating how given realizations behave 
for different values of the parameters. We present in 
Fig. [9] four different examples, for four different values of 
A, where four stochastic realizations (blue, cyan, purple, 
and pink lines) are compared with the classical trajec- 
tory (solid black line) . The top left panel corresponds to 
A — A0/2O. In this case, the noise is so small that the 
four stochastic realizations almost follow the classical tra- 
jectory. The top right panel corresponds to A = A^/IO, 
and the noise in the inflaton direction still plays no role in 
this case (i.e. the contour = 1 does not encompass the 
critical point). One sees that, when the trajectories enter 
the > 1 region, they start feeling the noise and that 
one of them (the pink realization) is even expelled to- 
wards the other global minimum. The bottom left panel 
corresponds to A = A^, and the noise is stronger as re- 
vealed by the "shaky" behavior of the realizations. The 
contour A^ = 1 appears but, clearly, the noise in the wa- 
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FIG. 10: Probability density functions in the {<j), tp) plane, at different times A^, starting from a Dirac distribution in the valley. 
The parameters chosen are A = 1.06347 Mp,, = M = 1.50398 Mp,, and fi = 7.74597 Mp,. 



terfall direction remains the main source of stochasticity. 
Finally, the bottom right panel corresponds to A = lOA^ 
and, this time, we are in a regime where the noise in the 
two directions is a priori important as indicated by the 
contours = 1 and = 1. This is especially clear for 
the cyan realization which climbs the inflationary valley. 
Therefore, it seems fair to say that our numerical code 
gives results that are completely compatible with elemen- 
tary expectations, which is an indication that it correctly 
calculates the behavior of the system. 

Finally, by simulating a high number of realizations, we 
have been able to calculate the correlation functions of 
various quantities of interest as well as their probability 



distributions. As an example, Fig.lTUlshows the probabil- 
ity density function in the (0, ip) plane, at different times 
iV, starting from a peaked distribution in the valley. One 
can see that the distribution is first roughly Gaussian 
and goes down the valley before setting over the critical 
point with "excrescences" growing towards the two min- 
imums of the potential, rendering the distribution highly 
non-Gaussian. Again, this plot confirms the previous dis- 
cussion and shows that the numerical codes used in this 
article are able to reproduce expected results in regimes 
where it is possible to guess (or to approximately calcu- 
late) the behavior of the system. 

In the following subsections, we present in more detail 
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FIG. 11: Distribution of the total number of e-folds realized 
during inflation. Classically, 505 e-folds are realized in the 
valley and 747 e-folds are realized during the waterfall stage, 
accounting for a total of 1252 e-folds. Clearly the mean value 
(A'^) ~ 50 is very different from the classical value which illus- 
trates well how important the stochastic effects in the vicinity 
of the critical point are. Despite this fact, it is also interest- 
ing to notice that, in the tail of the distribution, one can find 
realizations with a total number of e-folds larger than the 
classical value. 



our numerical results. 



D. Number of e- folds 

A first relevant well-defined physical quantity to study 
is the total number of e-folds realized during infiation 
since it provides a straightforward way to investigate the 
deviations from the classical picture. Of course, in order 
to calculate this quantity, one has to choose some initial 
conditions. Here, we take (pin/ipc — 1.0001 and ipin/M — 
'pj-^g parameters describing the shape of the po- 
tential are /.t/Mpi = 3190.4, M = cji^ = 0.1503Mp„ and 
A/Afp, = 0.01418. We see that this implies fiM/M^^ > 1 
and, therefore, we already know from the previous sec- 
tion that the number of e-folds during the waterfall phase 
will be large. As a matter of fact, it can be easily esti- 
mated upon using Eq. (|T7l). The above described choice 
is made in order to illustrate our point in the clearest 
way. It is important to stress that choosing other initial 
conditions would not drastically modify our conclusions. 
The classical prediction can be calculated in the slow-roll 
approximation using the formulas derived above, or using 
a numerical integration of the exact equations. It leads 
to a trajectory such that ^ 505 c-folds are realized in the 
valley and ^ 747 during the waterfall regime. The total 
number of e-folds is therefore ^ 1252. 

Then, we have computed the same quantity (for the 
same values of the initial conditions and of the param- 
eters) in the stochastic case. Obviously, for each real- 



FIG. 12: Distribution of the total number of e-folds realized 
during inflation, normalized by its maximum value, for differ- 
ent values of ii. The dependence on fi of {N) is consistent with 
the qualitative predictions of Sec. [Ill One also notices that the 
dispersion of the distribution increases with the number of e- 
folds realized in the flat region. This is due to the fact that 
the quantum effects broaden the distribution. 



ization one gets a different number, and the correspond- 
ing distribution is displayed in Fig. 1111 Let us now dis- 
cuss this figure. Probably, the most striking property of 
Fig.HHis that the distribution is peaked at a value which 
is completely different from the classical prediction. This 
clearly means that strong non-perturbative effects are at 
play. This also emphasizes the necessity of using a full 
numerical approach. Moreover, one sees that the stochas- 
tic contribution tends to diminish the total number of 
e-folds. This fact can be intuitively understood by notic- 
ing that most e-folds are realized in the region where the 
potential is very flat, around the critical point. Since this 
is precisely where the stochastic terms are dominant, the 
quantum kicks remove the system away from this region 
much faster than the classical roll; hence a lower number 
of c-folds is realized in this region. 

The tendency to escape faster from a region where the 
potential is very flat can be understood analytically on 
the example of small fleld inflation. In this single field 
model, the potential is given by 



V {ip) = 



1-1 



(62) 



where /x is a mass scale. Inflation proceeds from small to 
large values of the fleld. At the beginning of inflation, the 
potential is very flat and a large number of e-folds can 
be realized. For p — 2, the slow-roll and the perturbative 
Langevin equations can be integrated and solved exactly; 
see Ref. [4J]. Following Ref. |63j, one can then calculate 
the mean value of the total number of e-folds, 
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(63) 
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FIG. 13: Total number of e-folds (blue) and number of e- 
folds realized in the valley (green) as a function of A/Mn for 
the parameters indicated on the plot. The two numbers of e- 
folds are normalized to their classical counterparts. The solid 
lines represent the mean values of the distributions, while 
the colored surfaces represent the plus or minus 1 standard 
deviation areas. The vertical solid black line indicates the 
value of A^. 



In the present case and since </> ^ /i, one obtains 

4 / X 2 



{N) 



^ class 

In 



1 / M 
1927r2 \M^, 



(64) 



where A'ciass is the number of e-folds classically realized. 
From the above expression, it is clear that (TV) is smaller 
than iVciass since ip-m < ipend in this model. This result 
confirms the previous considerations: when the potential 
is very flat, the quantum kicks undergone by the inflaton 
field push it out of the flat region and, as a consequence, 
the total number of e-folds becomes smaller. 

Finally, it is also interesting to study how our results 
depend on the parameters of the model, especially on fi 
and A. As an example, Fig. [T^] shows the normalized 
distribution of the total number of e-folds for different 
values of the parameter /i. From Eqs. (fT6)) and (|47|) . one 
can see that the classical number of e-folds realized dur- 
ing both the inflationary valley and the waterfall regime 
increases with ^, which is consistent with the behavior of 
the mean values of the stochastic distributions observed 
in Fig. [T21 Moreover, the longer the field system stays 
in the flat region close to the critical point, the more its 
distribution gets stochastically broadened. This means 
that the dispersion of the distribution should evolve in a 
similar manner, which is exactly what is seen in Fig. [T^l 
Moreover, if one keeps increasing /z, one observes that 
the mean value of the distribution saturates at a value of 
~ 60 e-folds. 

As mentioned above, we have also studied how our 



FIG. 14: Number of e-folds realized between the exit of the 
stochastic regime (i.e. the moment where and A-,p are both 
smaller than 1) and the end of inflation. The cruxes stand for 
the mean values of the distributions, and the colored surfaces 
stand for the plus or minus 1 standard deviation areas. The 
horizontal dotted red line represents the total number of e- 
folds calculated in the absence of noise. The horizontal dotted 
orange line is the number of e-folds realized in the waterfall 
region in the absence of noise. Finally, the horizontal green 
dotted line represents the minimal number of e-folds required 
for inflation to be successful, i.e. 60. 



results depend on A. In the vacuum dominated regime, 
A is directly related to the energy scale of inflation. This 
quantity is constrained by the big bang nucleosynthesis 
and by the observations of the CMBR, namely. 



10^'^Mp, < A < 5 X lO^^Afp 



(65) 



The previous figures correspond to a regime where A > 
A^ (recall that, for the values of the parameters chosen 
here, we have A^ ~ 5.7 x 10~*Mpi) since we wanted to 
study a regime where the noise in the two-field directions 
is important. However, one can also wonder if the pre- 
vious conclusions, especially the fact that (N) < A^ciass, 
still hold in the regime A < A^. In Fig. [T31 we have 
computed the total number of e-folds as a function of A. 
The first thing we notice in this plot is that, for A > A^ 
(indicated by the vertical black line) , the total number of 
e-folds and the number of e-folds in the valley are smaller 
than their classical counterparts. This is of course com- 
patible with the previous considerations. In the regime 
10~^Afpi < A < A^, the number of e-folds in the valley is 
equal to its classical counterpart as expected since the in- 
flaton behaves classically but the total number of e-folds, 
hence the number of e-folds in the waterfall regime, is 



still smaller than N, 



So even in the absence of noise 



along the inflaton direction, the conclusion obtained be- 
fore remains valid. Finally, for A < 10~^Mp,, the noise 
is so small that the stochastic and classical number of 
e-folds are equal. 

In Fig. [TI] we have seen that, for A = 0.01418A/pi, one 
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FIG. 15: Level lines ei = 1 in the {(j), %p) plane. The pa- 
rameters chosen are the same as in the other plots. The red 
rhombus sits at the critical point, while the two green ones 
are located at the two minima. The definition of the an- 
gle 6 is shown for a trajectory leaving the inflationary region 
(i.e. crossing the ei = 1 line) at the black square. 



FIG. 16: Numerical distribution of the exit angle 9. The 
parameters are the same as in the other figures. The red line 
corresponds to the classical prediction. We notice that the 
distribution is extremely peaked, /S.6/9 ~ 10~^. 



has {N) ~ 50. Since the scales of astrophysical inter- 
est today left the Hubble radius during inflation about 
50 — 60 e-folds before the end of inflation, this would 
mean that we could have a direct observational window 
on the stochastic regime. In fact, this is not so because 
A = 1 means that H^/ei ~ 1. But H^/ei is precisely 
the overall normalization of the density perturbations 
power spectrum which is observed to be ~ 10"^. So, 
in fact, this shows that the value A — 0.01418Afp, is 
simply excluded by the CMBR measurements. For this 
reason, it is interesting to plot the number of e-folds per- 
formed between the moment the system becomes clas- 
sical (i.e. when and A^ are both smaller than 1) 
and the end of inflation. If this number is smaller than 
~ 50 — 60, this means that the corresponding value of 
A is excluded due to the above argument. The plot is 
represented in Fig. [Ml We see that for A > 10~^Mp,, 
the number mentioned above is indeed smaller than 60. 
All these values are therefore excluded. This means that 
for values of A such that A^ < A < lO^^Mpi, we are 
in a regime where the perturbations are not too large to 
be directly in contradiction with the CMBR (of course 
this does not guarantee that the correct normalization 
can be obtained) and where it is mandatory to take into 
account the stochastic effects in the two-field directions. 
If 10~^A/pi < A < A^, the stochastic effects dominate 
only in the tp direction, both in the valley and the wa- 
terfall phase. For lO'^Mpi < A < IQ-^Mp,, the water- 
fall regime becomes completely classical and, finally, for 
A < 10~^Mp,, the noise becomes so small that the full 
evolution in the valley and in the waterfall region can be 
described classically. 

Let us end this subsection with some remarks. We 
have seen that the stochastic number of e-folds is smaller 



than its classical counterpart as soon as A > lO^^Afp,. If 
A > A^ ~ 5.7 X 10~^Mp,, we are in a two-field regime. 
Moreover, if A > 10~^Mpi, the stochastic effects are so 
strong that the model is in contradiction with the am- 
plitude of the CMBR fluctuations. We expect these con- 
clusions to be very roughly independent of the choice of 
the other parameters, provided, of course, that one re- 
mains in the regime described in Sec. Ill Al In fact, to go 
further, one should explore the full parameter space, and 
one should carefully apply Cosmic Background Explorer 
(COBE) normalization to the model. When the waterfall 
regime plays an important role, this is not a trivial task. 

E. Inflation Exit Point 

Another relevant physical quantity is the inflation exit 
point, i.e. the location in the fleld space where inflation 
stops. The details of the subsequent (p)reheating phase 
strongly depend on these initial conditions, which are 
therefore important physical quantities [HI, [sl]. Infla- 
tion stops when the system crosses the ei = 1 level line 
in the field space. As a consequence, the exit point is 
necessarily located on this level line. It can be character- 
ized by the angle 9 between the line joining the closest 
minimum to the origin and the line joining this same 
minimum to the exit point. This parametrization and 
the £1 = 1 contours are represented in Fig. 1151 Since the 
stochastic effects taking place in the valley quickly render 
the distribution symmetrical in ip, the two minima are in 
fact put on an equal footing. Using the same method as 
before, one can calculate the classical prediction for this 
angle 9 and compare it with the corresponding stochastic 
distribution. The result is shown in Fig. 1161 

Here again, several comments are in order. First, un- 
like the distribution of the number of e-folds, the classi- 
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FIG. 17: Flow map of the classical slow-roll dynamics equa- 
tions. Regardless of their initial conditions, all the trajectories 
end up at the same exit point. 



cal prediction lies within the stochastic distribution. It is 
even more remarkable that the distribution is very nar- 
row. A priori, this is a surprising fact since the stochas- 
tic realizations in the vicinity of the critical point are 
extremely noisy, as can be seen in Fig. [HI Even if the 
stochastic trajectories are very different from one realiza- 
tion to another and spread over a large area in the field 
space, they eventually gather at the end of the waterfall 
phase to exit inflation at nearly the same point. We in- 
terpret this property in the following manner. Around 
the critical point, as was already mentioned before, the 
noise is quite strong and it quickly kicks the fields out 
of this region. As a consequence, the fields eventually 
land in a region where the noise is subdominant. There- 
fore, from that point, the fields will follow a classical 
trajectory. If there is a classical attractor, all the tra- 
jectories will converge towards this particular path, and 
inflation will always stop at the same point. This anal- 
ysis is confirmed by Fig. 1171 where we have plotted in 
the field plane the fiow lines of the classical equations of 
motion, Eqs. ([2]) and (jSj. As can be seen in the figure, 
after crossing the — 1 and = 1 level lines, which 
implies that the fields enter a region where the stochas- 
tic terms are subdominant, all the classical trajectories 
merge into a single one before crossing the ei = 1 level 
line and thus exiting the inflationary region. Then, one 
can check that this point corresponds to the angle 9 sin- 
gled out in Fig. [111 This classical attractiveness can also 
be formally established by studying the Lyapounov ex- 
ponent, in the direction orthogonal to the flow tangents. 
We conclude that, despite the strong quantum effects un- 
dergone by the flelds during the inflationary phase, the 
exit point is always the same (approximately, of course) 
in hybrid inflation. Moreover, this point turns out to be 
the classical one which provides a straightforward way to 
calculate its location. 



IV. CONCLUSION 



Let us now summarize our main findings. We have 
found that the quantum effects play an important role 
in hybrid infiation, especially in the vicinity of the criti- 
cal point (this seems to be a general feature of multiple 
field models of inflation, as soon as flat directions are 
present in the potential). As a consequence, the classical 
picture presented in Sec.|lT]has to be substantially modi- 
fied. This can be done in the framework of stochastic in- 
fiation, where the inflationary dynamics is driven by two 
coupled Langevin equations. Given that the stochastic 
effects can be strong in the two directions in fleld space, 
we have used a numerical approach to solve these equa- 
tions. Then, we have derived the distributions of two 
relevant quantities, namely the total number of e-folds 
realized during inflation and the exit point. We have 
shown that, when the stochastic noise plays a role, the 
distribution of the number of e-folds is peaked at a value 
which is different from the classical prediction. This is 
due to the fact that, in the neighborhood of the critical 
point, the potential is very flat and the quantum kicks 
quickly move the system away from this region. On the 
other hand, the distribution of the exit point of inflation 
leads to conclusions which are apparently at odds with 
this picture since it is extremely peaked over the classi- 
cal prediction. But, in fact, this property is due to the 
attractiveness of the classical flow and is not at all in 
contradiction with the previous considerations. 



An important question that remains to be addressed 
in more detail is the impact of our results on the ob- 
servable predictions of hybrid inflation. For instance, it 
would be of utmost importance to study how the quan- 
tum effects can modify the power spectra of cosmological 
fluctuations. It was recently emphasized in Ref. [2^ that 
hybrid inflation can lead to a red spectrum, and it would 
be interesting to investigate the influence of the quan- 
tum effects on t his p rediction, both for the adiabatic and 
entropy modes 1665. Also, the explicit computation of 
the probability density functions in the fleld space pro- 
vides us with a means to calculate the non-Gaussianities 
of this model. We intend to come back to those issues in 
the future. Maybe the most important conclusion of our 
work is that the richness of multiple inflation - namely, 
the presence of entropy modes a priori produced during 
the waterfall regime, the highly non-trivial phase of pre- 
heating, the strong quantum effects - implies that it is 
not simple to derive the corresponding observable pre- 
dictions and that, most of the time, these ones cannot 
be obtained in a simple single field effective model. This 
is an important conclusion that one should keep in mind 
when analyzing the future high accuracy CMB data and 
their implications for inflation. 



16 



Acknowledgments 

We would like to thank S. Clesse and C. Ringeval for 
useful discussions. 



A. A. Starobinsky, Phys. Lett. B91, 99 (1980). [34 
A. H. Guth, Phys. Rev. D23, 347 (1981). [35 
A. D. Linde, Phys. Lett. B108, 389 (1982). 
A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, [36 
1220 (1982). [37 
A. D. Lindo, Phys. Lett. B129, 177 (1983). 
J. Martin, Braz. J. Phys. 34, 1307 (2004), astro- [38 
ph/0312492. 

J. Martin, Lect. Notes Phys. 669, 199 (2005), hep- [39 
th/0406011. 

J. Martin, Lect. Notes Phys. 738, 193 (2008), 0704.3540. [40 
V. F. Mukhanov and G. Chibisov, JETP Lett. 33, 532 [41 

(1981) . [42 
V. F. Mukhanov and G. Chibisov, Sov. Phys. JETP 56, 
258 (1982). [43 
S. Hawking, Phys. Lett. B115, 295 (1982), revised ver- 
sion. [44 
A. A. Starobinsky, Phys. Lett. B117, 175 (1982). 
A. H. Guth and S. Y. Pi, Phys. Rev. Lett. 49, 1110 [45 

(1982) . 

J. M. Bardeen, P. J. Steinhardt, and M. S. Turner, Phys. [46 
Rev. D28, 679 (1983). 
J. Martin and C. Ringeval, JCAP 0608, 009 (2006), [47 
astro-ph/0605367. 
L. Lorcnz, J. Martin, and C. Ringeval, JCAP 0804, 001 [48 
(2008), 0709.3758. 

L. Lorcnz, J. Martin, and C. Ringeval, Phys. Rev. D78, [49 
063543 (2008), 0807.2414. 
J. Martin and C. Ringeval, Phys. Rev. D82, 023511 [50 

(2010), 1004.5525. 

J. Martin, C. Ringeval, and R. Trotta, Phys. Rev. D83, [51 
063524 (2011), 1009.4157. 

A. D. Linde, Phys. Rev. D49, 748 (1994), astro- [52 
ph/9307002. [53 
E. J. Copeland, A. R. Liddle, D. H. Lyth, E. D. Stew- [54 
art, and D. Wands, Phys. Rev. D49, 6410 (1994), astro- 
ph/9401011. ' [55 

E. Halyo, Phys. Lett. B387, 43 (1996), hep-ph/9606423. 
P. Binetruy and G. DvaU, Phys. Lett. B388, 241 (1996), [56 
hep-ph/9606342. [57 
G. Dvali, Q. Shafi, and R. K. Schaefer, Phys. Rev. Lett. 
73, 1886 (1904), hcp-ph/9406319. [58 
R. Kallosh and A. D. Linde, JCAP 0310, 008 (2003), 
hep-th/0306058. [59 
S. Clesse and J. Rocher, Phys. Rev. D79, 103507 (2009), 
0809.4355. 

S. Clesse, Phys. Rev. D83, 063518 (2011), 1006.4522. [60 
A. A. Abolhasani, H. Firouzjahi, and M. H. Namjoo, 
Class. Quant. Grav. 28, 075009 (2011), 1010.6292. 
A. Vilenkin, Nucl. Phys. B226, 527 (1983). [61 
A. A. Starobinsky (1986). 

A. Goncharov, A. D. Linde, and V. F. Mukhanov, Int. J. [62 
Mod. Phys. A2, 561 (1987). 
Y. Nambu and M. Sasaki, Phys. Lett. B205, 441 (1988). [63 
Y. Nambu and M. Sasaki, Phys. Lett. B219, 240 (1989). 



H. E. Kandrup, Phys. Rev. D39, 2245 (1989). 

K.-i. Nakao, Y. Nambu, and M. Sasaki, Prog. Theor. 

Phys. 80, 1041 (1988). 

Y. Nambu, Prog. Theor. Phys. 81, 1037 (1989). 

S. Mollerach, S. Matarrese, A. Ortolan, and F. Lucchin, 

Phys. Rev. D44, 1670 (1991). 

A. D. Linde, D. A. Linde, and A. Mezhlumian, Phys. 

Rev. D49, 1783 (1994), gr-qc/9306035. 

A. A. Starobinsky and J. Yokoyama, Phys. Rev. D50, 

6357 (1994), astro-ph/9407016. 

A. D. Linde, Phys. Lett. B175, 395 (1986). 

A. D. Linde, Mod. Phys. Lett. Al, 81 (1986). 

J. Garcia-BoUido, A. D. Linde, and D. Wands, Phys. Rev. 

D54, 6040 (1996), astro-ph/9605094. 

J. Martin and M. Musso, Phys. Rev. D73, 043516 (2006), 

hep-th/0511214. 

J. Martin and M. Musso, Phys. Rev. D73, 043517 (2006), 
hep-th/0511292. 

A. M. Green and A. Mazumdar, Phys. Rev. D65, 105022 
(2002), hep-ph/0201209. 

E. D. Stewart and D. H. Lyth, Phys. Lett. B302, 171 
(1993), gr-qc/9302019. 

A. R. Liddle, P. Parsons, and J. D. Barrow, Phys. Rev. 
D50, 7222 (1994), astro-ph/9408015. 
S. Groot Nibbelink and B. van Tent, Class. Quant. Grav. 
19, 613 (2002), hop-ph/0107272. 

R. Easthor and J. T. Giblin, Phys.Rev. D72, 103505 
(2005), astro-ph/0505033. 

D. J. Schwarz, C. A. Terrero-Escalante, and A. A. Garcia, 

Phys. Lett. B517, 243 (2001), astro-ph/0106020. 

S. M. Leach, A. R. Liddle, J. Martin, and D. J. Schwarz, 

Phys. Rev. D66, 023515 (2002), astro-ph/0202094. 

S. Clesse, AIP Conf. Proc. 1241, 543 (2010), 0910.3819. 

S. Clesse (2010), 1006.4435. 

S. Clesse, C. Ringeval, and J. Rocher, Phys. Rev. D80, 
123534 (2009), 0909.0402. 

A. A. Abolhasani, H. Firouzjahi, and M. Sasaki (2011), 
1106.6315. 

D. H. Lyth (2010), 1005.2461. 

A. A. Abolhasani and H. Firouzjahi, Phys. Rev. D83, 
063513 (2011), 1005.2934. 

H. Kodama, K. Kohri, and K. Nakayama, 
Prog.Thcor.Phys. 126, 331 (2011), 1102.5612. 

F. FincUi, G. Marozzi, A. A. Starobinsky, G. P. 
Vacca, and G. Venturi, Phys. Rev. D79, 044007 (2009), 
0808.1786. 

F. Finelh, G. Marozzi, A. A. Starobinsky, G. P. 
Vacca, and G. Venturi, Phys. Rev. D82, 064020 (2010), 
1003.1327. 

F. Finelli, G. Marozzi, A. Starobinsky, G. Vacca, and 

G. Venturi (2011), 1102.0216. 

K. Enqvist, D. G. Figueroa, and G. Rigopoulos (2011), 
1109.3024. 

S. Gratton and N. Turok, Phys. Rev. D72, 043507 
(2005), hep-th/0503063. 



17 



[64] J. Garcia-Bellido and A. D. Linde, Phys. Rev. D57, 6075 [66] L. P. Lcvassour, G. Laportc, and R. Brandenberger, 
(1998), hcp-ph/9711360. Phys. Rev. D82, 123524 (2010), 1004.1425. 

[65] F. FincUi and R. H. Brandenberger, Phys. Rev. D62, 
083502 (2000), hep-ph/0003172. 



